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Abstract: It is argued that experimental constraints on theories of asymmetric dark matter 

►^ (ADM) almost certainly require that the DM be part of a richer hidden sector of interacting 

'^ states of comparable mass or lighter. A general requisite of models of ADM is that the vast 

QQ majority of the symmetric component of the DM number density must be removed in order 

'^ to explain the observed relationship O^b ~ ^^dm via the DM asymmetry. Demanding the 

C^ efficient annihilation of the symmetric component leads to a tension with experimental limits 

^vq if the annihilation is directly to Standard Model (SM) degrees of freedom. A comprehensive 

'""' effective operator analysis of the model independent constraints on ADM from direct detection 

^ experiments and LHC monojet searches is presented. Notably, the limits obtained essentially 

S^ exclude models of ADM with mass 1 GeV < ttidm ^ 100 GeV annihilating to SM quarks 

H via heavy mediator states. This motivates the study of portal interactions between the dark 

and SM sectors mediated by light states. Resonances and threshold effects involving the new 

light states are shown to be important for determining the exclusion limits. 
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1 Introduction 

Asymmetric dark matter (ADM) provides a compelling framework for explaining the existence 
and comparable magnitudes of the present day baryon density i^B and the dark matter (DM) 
relic density ^bm- In traditional cosmology the origins of Qb and f^DM are distinct and 
unrelated. Baryogenesis is conventionally attributed to CP violating processes in baryons 
or leptons, whilst J^dm is assumed to be determined by particle species decoupling from 
the thermal bath (freeze-out) [1-3], or out-of-equilibrium species moving towards thermal 
equilibrium with the visible sector (freeze-in) [4]. Thus, a priori, the DM and baryon densities 
could have differed by orders of magnitude and the notably close proximity f^oM ~ 4.56 x Qb 
has motivated a concentrated research effort to connect the genesis mechanisms for baryons 
and DM [5-16]. 

For the DM to be stable and linked to the baryon asymmetry, in addition to the ob- 
served quantum number which stabilises the state (e.g. Z2 i?-parity), the DM must carry an 
additional conserved quantum number analogous to baryon number, which we denote X. At 
low energies U(1)b and U(l)x are approximate global symmetries of the visible and hidden 
sector, respectively. However, at higher energies these are no longer good symmetries and the 
only conserved quantity is B — L + X. Consequently, models can be constructed in which 
asymmetries in baryon, lepton and DM number are linked, and thus provide a natural ex- 
planation for the relative coincidence of the cosmic relic densities. There now exist a large 
variety of mechanisms for connecting DM with baryogensis [17-38], but they roughly fall 
into two broad categories: sharing and cogenesis. In models of sharing one begins with a 
pre-existing asymmetry in some (approximately) conserved quantum number and postulates 
a connector operator through which the asymmetry may be transferred among B — L and 



X. In the oogenesis scenario one assumes that there is no initial asymmetry and that the 
asymmetries in baryons and DM are simultaneously generated via some out-of-equilibrium 
CP violating process [34-38] . In many models of oogenesis the baryon and DM asymmetries 
are equal and opposite, such that there is no overall asymmetry in the quantity B — L + X. 
Since natural realisations of ADM lead to comparable asymmetries in B and X, to account 
for Qb/^dm ~ 5, the DM mass ttidm must be of similar magnitude to the proton mass 
nip ss 0.94 GeV. We shall refer to the range 1 GeV < tbdm ^ 10 GeV as the ADM mass 
range and focus much of our attention on this mass region. Although the ADM mass range 
is preferred theoretically, we include limits on model of ADM with DM mass in the range to 
1 GeV < moM < 10^ GeV. 

In standard theories of baryogenesis an asymmetry r]B is generated between the baryons 
and anti-baryons which ultimately leads to the visible universe. Due to the strong nuclear 
force the symmetric component of the baryon density pair-wise annihilates, leaving only 
the asymmetric component of the baryons, resulting in the matter dominated universe we 
observe. The framework of ADM proposes that the (non-self-conjugate) DM undergoes a 
similar cosmological history to the baryons and the DM relic density is due to an asymmetry 
between the DM and its antiparticle. Similar to baryons, the DM must have sufficiently 
large interactions such that only a negligible symmetric component remains, in which case 
the baryon and DM relic densities are related by 

^DM ^ ^B- (1-1) 

rripTjE 

To link Qb and 0dm in this manner it is essential that the DM relic density be due to the 
DM asymmetry and consequently, as discussed in [35], it is a generic requirement that the 
vast majority of the symmetric component is removed. 

In this paper we shall remain agnostic as to the source of the DM asymmetry, instead 
focusing on a general feature of ADM from which model independent bounds can be derived. 
Whilst DM annihilation is more efficient for larger values of the inter-sector couplings, it is 
precisely these inter-sector couplings which are constrained through direct detection and 
collider searches. Hence, the need for efficient annihilation results in a tension between 
cosmological requirements and direct search constraints. Generically, the symmetric part 
of the DM will be several orders of magnitude larger than the asymmetric component, as 
the ratio of asymmetric to symmetric DM is of comparable magnitude to the CP violating 
parameter determining the asymmetry and is typically < 10~^. As a result, these constraints 
are sufficiently strong to rule out large regions of parameter space. 

Since momentum transfer is low in direct detection experiments, mediators with masses 
~ 100 MeV or greater are sufficiently heavy to be integrated out. Consequently the variety of 
portal interactions between the hidden and visible sectors may be parameterised via effective 
contact operators. We shall demonstrate that the current experimental limits exclude ADM 
models in which the DM couples to SM quarks with minimal flavour structure via a heavy 
mediator, for 1 GeV < /tt-dm ^ 100 GeV . However, effective operators do not provide a good 
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Table 1. Set of effective operators O of dimension 6 or less connecting SM quarks / to scalar (f> or 
fermion ^ DM with universal couplings. We consider 3 additional natural contact operators O with 
quark mass dependent couplings, note O = ^O. It is indicated whether the corresponding direct 
detection cross-section is DM velocity v or momentum transfer q suppressed and if the interaction 
couples to nuclei in a spin-dependent (SD) or independent (SI) manner [39]. 



description of the collider limits and the relic density calculation for mediators with masses 
< 100 GeV and, hence, in Section 4 we examine the effects of such light mediators. In Table 1 
we list all effective operators up to dimension 6 connecting SM quarks / to complex scalar (p or 
Dirac fermion t/j DM.^ For definitiveness and to avoid stringent flavour-changing constraints 
we will typically assume that the contact operators couple universally to all quark flavours. 
In the case of scalar and pseudoscalar interactions we shall also study the case where the 
effective coupling to SM quarks is proportional to the relevant quark mass, which is arguably 



more natural, and we denote these m„-dependent operators 



'-0. Operators 0$, O 



^Since the DM must carry a conserved \]{l)x Majorana fermions or real scalars are not suitable candidates 

e'^'''"'apcr, the operators tpia'^'^'y^tpfa^^f and 



t^^^A 



for ADM. Recall a^"' = ^{j^Y - Yl*") and since cr'"''7 

ijia'^^'ipfaf^vf are equivalent to ipia^'' ^jj f a fj^v'y^ f and tpo-^'''y'^'ipfaf_iv'y''f, respectively. 



Ops and Ofp correspond to CP violating interactions and one might expect their coefficients to 
be suppressed relative to those of the other, CP conserving, operators, consequently we shall 
typically omit these contact interactions for brevity. We have checked that the requirements 
for efficient annihilation of the symmetric component and experimental limits for these CP 
violating operators are comparable to examples which will be presented in detail. 

The direct detection cross-section of certain operators are suppressed by the small DM 
velocity v ~ 10~^c or momentum transfer \q\ = \pi — pf\ < 0.1 GeV. For operators with 
suppressed scattering cross-sections the current direct detection limits are insufficient to pro- 
vide any useful constraints on the DM over the whole mass range. However, the high energy 
production cross-sections remain unsuppressed and, consequently, monojet searches provide 
the leading limits on such operators. Approaches based on effective operators have been em- 
ployed previously to study the constraints from direct detection and colliders on conventional 
DM models [39-51], however, the existence of an asymmetry in the hidden sector leads to 
deviations in the relic density calculations [52]. Whilst there have been previous papers on 
contact operators in the context of ADM [53, 54], in the case of the former the effects of the 
DM asymmetry are not fully accounted for, while the latter work considers only a limited 
selection of possible operators. In the present work we consider the full set of contact oper- 
ators (correcting some numerical errors contained in previous works) and include the latest 
experimental constraints from both direct detection and LHC monojet searches (including 
the recent CMS results with 4.67fb~ ). Moreover, in later sections, we go beyond the contact 
operator analysis and demonstrate that if the DM is connected to SM quarks via exchange 
of a light mediator state then resonance and threshold effects become very important in 
determining the exclusion limits. 

As discussed previously, if the DM annihilates directly to SM quarks then there is a 
tension between the efficient removal of the symmetric component and experimental searches. 
However, if the DM is embedded into a richer hidden sector, with a spectrum of lighter 
states into which the symmetric component can annihilate (which later decay to the visible 
sector), then this severs the connection between direct search constraints and cosmological 
requirements. As we show that effective theories of ADM with ttt-dm ~ "ip and minimal 
flavour structure which decay directly to SM quarks are excluded for all types of connector 
operator, this leads to the striking conclusion that in order to realise ADM naturally, the DM 
must generally be accompanied by additional hidden states of comparable mass or lighter, 
either to mediate the interactions to the visible sector, or into which the symmetric component 
of the DM may annihilate. The fact that ADM must be part of a larger hidden sector is an 
important realisation, since this will generally lead to significant changes to the cosmology 
and phenomenology, which we shall address in a companion paper [55]. 

This paper is structured as follows, in Section 2 we review the standard treatment for 
calculating the DM relic density in the presence of an asymmetry. In Section 3 we present a 
comprehensive study of effective operators and compare the requirements for successful ADM 
models to the experimental limits from direct detection and monojet searches. The results of 
Section 3 strongly motivate a careful study of light mediators and we present an analysis of 



the scalar and pseudoscalar portal interactions due to light mediators in Section 4. Finally, 
in the concluding remarks we discuss the implications these results have on model building 
with ADM and emphasise that the null experimental searches strongly imply that successful 
models must feature an extended hidden sector. 



2 Relic abundance with an asymmetry 

The DM asymmetry alters the relic density calculation from the conventional case without 
an asymmetry. Relic density calculations in the ADM paradigm were recently discussed in 
[56, 57] and we shall follow these analyses in determining the relic density of the symmetric 
component. Generically, it was found in these analyses that the annihilation cross-section is 
required to be larger than in the case with no asymmetry. 

We define the standard variable x = m/T, for a particle of mass m at temperature T 
and denote by xp the inverse scaled decoupling temperature of DM and antiDM, assuming 
that the only X number changing interactions are pair-annihilations of the DM and its anti- 
partner. Following [56], the yields for X and X in the limit that re — )• oo are given by: 

f/DM 1^ ??DM 



^DM 



Y, 



where t/dm 



1 - exp [-r/DM-''w] ' ' ^^ exp [qnuJ^] - 1 ' 

^DM — ^%M ^^ ^^^ asymmetric yield. The quantity uj is defined by 
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mnuMpiy/g^ 



(2.1) 



(2.2) 



where w-dm is the DM mass, Mpi = 2.4 x 10^^ is the reduced Plank mass and g^:{T) is the 
number of available relativistic degrees of freedom at a given temperature T. The annihilation 

integrals J are given by 

{av) 



J 



dx. 



Xp 



Making an expansion of the annihilation cross-section" 
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allows us to express eqn. (2.1) as follows 
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Using the tables from DarkSUSY [58] for the temperature dependence of g*{T), we model 
the behaviour of xp via the approximation introduced in [56] : 
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(2.6) 



As discussed in Section 4 and the Appendix, this expansion is not appropriate in the case of hght mediators 
and must be replaced with an alternative analysis and numerical methods. 



where xpo is the inverse scaled decouphng temperature in standard freeze-out with no asym- 
metry. The present rehc density is given by 



nBuh^ = 2.76 X 10^ (??: 
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where Isym = 21j)5j is the yield of the symmetric component and h = 0.710 it 0.025 is 
the Hubble constant in units of 100 kms^^Mpc^ . In the limit ?7dm — ^ 0, one recovers 
the standard symmetric DM expressions. The observed value of the DM abundance [59] is 
^DM^^ = 0.1109 lb 0.0056. In order to realise the requirement that the asymmetric com- 
ponent constitutes the majority of the present relic density, we demand that the symmetric 
component composes < 1% of the relic density, via the constraint: 
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(2.8) 



Varying the constraint on the symmetric component in a wide range around 1% does not 
materially affect our conclusions concerning the required annihilation cross-section. 

The annihilation cross-sections for the operators in Table 1 connecting fermion DM and 
SM quarks to 0{v'^) are ^ 
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^These expressions correct previous errors in the literature and have been checked both by hand and with 
Feyncalc [60]. 
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while the annihilation cross-sections for operators with scalar DM are 
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In the above, the cross-sections corresponding to scalar or pseudoscalar interactions {Of , Op 
and of) have been provided for the case of universal coupling, the cross-sections for the niq- 
dependent couplings may be obtained by multiplying the cross-section (under the summation) 
by a factor of rni/ h? . To illustrate that the omitted CP violating operators are comparable 
to the CP conserving interactions we have included the operator Oz,^. 

These expressions are used to derive the maximum value of the scale A permitted to 
ensure a suitable depletion of the symmetric component for each operator as a function of 
"1-DM) where we assume that only one operator is turned on at a time. Imposing the constraint 
on the symmetric abundance, eqn. (2.8), the results obtained are presented graphically in 
Figures 1-3. The black curves present upper bounds for A for each operator as a function of 
fnuM- The resulting relic density curves may then be constrained through direct detection 
and collider limits, which we derive in the next section. 



3 Direct search constraints on contact operators 

The constraints on each contact operator in Table 1 depend strongly on whether the cor- 
responding scattering cross-section is spin- independent or spin-dependent. If the scattering 
cross-section is suppressed by powers of v or q then direct detection limits are severely weak- 
ened and, in fact, in all cases do not present any significant constraint on the allowed DM 
mass range. Rather, for the suppressed operators the leading bounds come from collider 
monojet searches, which we shall discuss shortly. 

The scattering cross-section per nucleon for each of the non-suppressed spin-independent 
interactions under consideration which couple universally to SM quarks^ are [40, 41] 



of _ _•_ 2f2 ot ^ JL 2 



(3.1) 



where u„ = / ™°^V^ n is the nucleon-DM reduced mass and L is the DM effective couplings 
to protons. For consistency with later calculations, we use the default values for the factor fp 
from micrOMEGAs [61]. The non-suppressed scattering cross-sections for the spin-dependent 
operators of interest with universal couplings are given by 

4l^^^-sl^ih^^i\y>l\ ■ (3-2) 

The Ag account for the spin content of the nucleon and we use the values derived in [62] from 

\ 2 
,p> 




the COMPASS [63] results: (Eg^g) ~ 0.32. Note that in the non-relativistic limit the 
form of the tensor operators qa^^q coincides with the axial-vector operator q^^^^^q up to a 
numerical factor [43]. From the scattering cross-sections given above we can use the current 
direct detection limits to place a lower bound on the scale A. 

The lower bounds placed on A by direct detection are displayed in Figures 1-3 as solid 
coloured curves. The leading limits on spin-independent interactions from direct detection 
experiments for rriDM ^ 1 GeV come from CRESST [64] (orange), while in the intermediate 
mass range 1 GeV< m-oM ^10 GeV the dominant constraints come from a combination of 
XenonlO [65] (green) CDMS [66] (cyan) and DAMIC [67] (magenta). For moM > 10 GeV 
the leading spin- independent limits are provided by the XenonlOO experiment [68] (blue). 
The strongest direct detection constraints on spin-dependent interactions are from Simple 
[69] (Stage 2: light purple; Combined: dark purple) for rriDM > 5 GeV and from CRESST 
[64] for hghter DM. 

A complementary set of constraints on interactions between SM quarks and DM comes 
from collider experiments, in particular the recent LHC searches for jets and missing energy 



*SI cross-section of operators which couple proportional to niq may be obtained by suitable rescalings; 
where we have assumed that fn ~ fp- 
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Figure 1. Limits on the scale A for the effective operators from direct detection and monojets for 
scalar DM (f) and fermion DM ij), as a function of DM mass mDM- The black curve corresponds to the 
minimum annihilation cross-section necessary to reduce the symmetric component to 1%. Constraints 
are from XenonlOO (blue), XenonlO (green), CDMS (cyan), CRESST (orange), DAMIC (magenta) 
and ATLAS lfb~^ (red, dashed) and CMS 4.67 fb"""^ (blue, dashed) monojet searches. For ttidm ^ A 
effective operators no longer provide a good description and can not be reliably used to calculate the 
relic density requirements. Moreover, the monojet limits are no longer reliable much below A < 100 
GeV due to the experimental cuts employed by ATLAS and CMS. Viable models of ADM employing 
the listed effective operators must lie in the shaded parameter regions. Note that contact operators due 
to scalar mediators are studied for both universal couplings to quarks and rriq-dependent couplings. We 
see that, except for the operator j^(j)^(j)qq around ttidm ~ 1 GeV, successful models of ADM involving 
contact operators are excluded for 1 GeV < todm ^ 100 GeV, which includes the range in which ADM 
is most well motivated. 
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Figure 2. Limits on A for operators witlr spin-dependent direct detection cross-sections, viable 
parameter regions are shaded. Constraints are from Simple (Stage 2: light purple; Combined: dark 
purple), CRESST (orange), ATLAS Ifb"^ (red, dashed) and CMS 4.67 fb"^ (blue, dashed). 






Figure 3. Limits on A for operators with v ox q suppressed direct detection cross-sections, viable 
parameter regions are shaded. Limits are from ATLAS 1 fb^ (red, dashed) and CMS 4.67 fb^ (blue, 
dashed). The interesting ADM range ttidm ^ 10 GeV is excluded in all cases and, with the exception 
of the ^tp'j^ipq'j^q operator, this exclusion extends up to moM ^ 100 GeV. 



10 



[70, 71].^ Operating at ^/s = 7 TeV, the ATLAS and CMS collaborations have reported no 
apparent excess in the production of jets and missing energy with an integrated luminosity 
of Ifb^ and 4.67fb~ , respectively. The high px ATLAS analysis [70] selected events with 
a primary jet with px > 250 GeV and pseudorapidity \rj\ < 2, and $x > 220 GeV, events 
with a secondary jet with px < 60 GeV and |?7| < 4.5 were also permitted. The collaboration 
reported 965 events. The SM prediction for this process is 1010 it 75 events. The ATLAS null 
result excludes any new contributions to the production cross-section greater than 0.01 pb. 
The ATLAS monojet searches have been previously used to place limits on certain models of 
standard symmetric DM [47-51]. Here we use the LHC monojet bounds to derive constraints 
on contact interactions between SM and DM states for the ADM scenario and employ the 
most recent results. The recent CMS search [71] excludes new contributions to the production 
cross-section greater than 0.02 pb and presents a comparable set of constraints to the ATLAS 
1 fb~ limits. The CMS analysis considered events with a primary jet with px > 110 GeV 
and pseudorapidity \7]\ < 2.4, and $x > 350 GeV, they also allowed events with a secondary 
jet with Px > 30 GeV, provided the azimuth angle difference of the jets satisfied Aip < 2.5. 
CMS observed 1142 events in good agreement with the SM prediction of 1224 it 101 events. 

To calculate the limits on contact operators coupling DM to quarks, we use the program 
CalcHEP [72]. The SM Lagrangian is supplemented with each operator in separate instances 
and we study the total cross-section pp — )• jXX where X is the DM state and j is a jet. To 
model the ATLAS search we apply the following cuts on the events ^x > 220 GeV, px > 250 
GeV, and pseudorapidity |??| < 2 and for CMS we take $x > 350 GeV, px > HO GeV, and 
pseudorapidity |?7| < 2.4. We do not consider events with additional jets in either cases. With 
these cuts we may reliably assume that the efficiency of the searches is close to 100%. For 
each operator we determine the minimum value of A for which the new contribution to the 
production cross-section is permitted by monojet searches. Since the contact operators are 
generally suppressed by multiple powers of A, 0(1) changes to the cross-section result in only 
small deviations in the limit on the scale A. The resulting limits are plotted in Figures 1-3 
as dashed blue (CMS) and red (ATLAS) curves and supplement the constraints from direct 
detection. 

For rriBM ^ A the contact operators no longer give a good description and one is required 
to consider the UV completion of the effective theory. Moreover, because of the cuts applied to 
the events by the experimental collaborations, the effective operators do not provide reliable 
bounds for monojet searches if the limit on A is much below the px cut ~ 100 GeV. This is the 
case with the operator ^(j)(l>^qq, and the fermion DM, scalar and pseudoscalar ?TT,g-dependent 
operator monojet bounds are close to the limit at which the effective theory breaks down. 
Consequently, these portal interactions should really be studied in the light mediator regime, 
which we discuss in detail in Section 4. Note, however, that the Tevatron monojet searches 
utilise much lower momentum cuts, requiring a primary jet with px > 10 GeV, \rj\ < 1.1, 
and $x > 60 GeV (90% efficiency), and thus the contact operator description for monojet 



^Note that LHC monophoton searches provide sUghtly weaker bounds to those from monojets [71]. 
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searches are reliable for A > 10 GeV. Moreover, results from CDF [73] lead to constraints 
which are roughly comparable to the CMS limits for 1 GeV < w-dm ^ 10 GeV. 

To correctly interpret Figures 1-3, recall that the black curves provide an upper bound on 
A to ensure suitably efficient annihilation of the symmetric component, whilst the coloured 
curves present lower bounds from searches. A given operator may present a viable ADM 
model only if there exists some mass region for which a sufficiently low A has not been 
excluded by direct search constraints. It is immediately apparent that there exists a tension 
between the experimental constraints and the upper bound on A required for the efficient 
annihilation of the symmetric component over much of the model space. 

ADM most naturally resides in the mass range 1 GeV < m^M < 10 GeV. From the 
plots in Figures 1-3 we see that essentially no operators are viable in this preferred mass 
range. For scalar DM a small region of parameter space remains at tjidm ~ 1 GeV for the 
operator -^4>(j)^qq. All of the other spin-independent non-suppressed operators are excluded 
by direct detection experiments up to at least ttidm = 100 GeV. The strongest bounds on the 
remaining fermion DM operators in the ADM mass region are from collider searches, which 
are sufficient to exclude these operators over the mass range 1 GeV < tjidm ^ 10 GeV. With 
a small number of exceptions the common bound on the scale A due to monojet searches is 
roughly A > 1 TeV for muM ^ 100 GeV and present stronger bounds than direct detection 
constraints in several cases. For the non-suppressed spin-independent operators, since the 
ADM mass range is already excluded by direct detection experiments, these collider limits 
do not impose significant new constraints. In contradistinction the LHC searches provide 
virtually the sole limits on spin-dependent interactions and the velocity suppressed operators 
over the mass region of interest. 

On the other hand, inspecting Figures 1-3 we note that TeV scale ADM (in distinction 
to ADM in the natural mass range 1 GeV < w-dm ^ 10 GeV) allows most effective operators 
to remove the symmetric component efficiently without conflict with current detection limits. 
We have previously argued that if the asymmetries are equal, then DM should be roughly 5 
GeV, however the TeV scale provides an alternative mass range for which the correct value 
of f^DM can he realised [5, 8]. This is because for large masses DM production is Boltzmann 
suppressed by a factor exp (— '"j^^ ) for todm > T. But, the TeV scale ADM scenario is 
rather less appealing, since the DM relic density is exponentially sensitive to changes in 
rriDM- Alternatively, in sharing models the operator which transfers the asymmetry may be 
inefficient, and this can be used to generate large discrepancies between the baryon and DM 
asymmetries. Whilst such constructions permit a wider range of m-oM, it as the mass of the 
DM rises it becomes increasingly difficult to explain the coincidence JIdm/^b ~ 5 via the DM 
asymmetry. Thus, although a wider range of ADM masses are certainly worth contemplating, 
since a natural expectation in many models of ADM is that the asymmetries in X and B 
should be comparable, we have focused much of our discussion on the case ttt-dm ~ ""T-p- 

A further source of limits on contact operators come from invisible quarkonium decays. 
Clearly, only bound states heavier than 2mDM can decay invisibly to the DM states. The 
principle constraints on the ADM region come from decays of the T(15) mesons, as measured 
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at the B-factories BaBar and CLEO [74]. However, we find that the Umits on the contact 
operators of Table 1 coming from these invisible decay searches are not competitive with pre- 
vious bounds from monojet searches or direct detection. There are also additional constraints 
on ADM from astrophysics [75-79]. In particular for ADM models with scalar DM there are 
additional constraints from old compact stars [80-82]. These limits are sufficiently strong 
to exclude scalar ADM over a large mass range (including the natural ADM mass range), 
however this conclusion can be circumvented with further model building. For instance, these 
limits no longer apply if the scalar DM is not fundamental, but a composite state of fermions 
due to some new hidden sector strong dynamics, such that Fermi repulsion of the constituent 
states becomes important before the Chandrasekhar limit is reached [82, 83]. 

In summary, the combination of direct detection and LHC limits excludes DM-quark 
contact interaction in the ADM mass region for all contact operators with minimal flavour 
structure. Moreover, with the exception of the pseudoscalar operator ^tp'j^fpq'y^q, all of these 
contact operators are disallowed by experimental searches for 1 GeV < wdm ^ 100 GeV. 
The conclusion is striking: models of ADM with the DM at the natural mass scale (1-10 GeV) 
decaying directly to SM quarks via heavy mediators is completely excluded and only light 
mediators provide suitable portal interactions. The case of light mediators is not a trivial 
extension of the effective operator analysis since the results for the monojet and relic density 
calculations change drastically. We shall discuss the limits and model building opportunities 
of light mediators in the next section. 

4 Light mediator analysis 

A proper treatment of light mediators requires a careful analysis which includes effects of 
resonances and mass thresholds. Mediators with mass comparable to, or less than, the mo- 
mentum cuts employed in the collider searches can lead to a substantial weakening of the 
monojet bounds. Also, as noted previously, contact operators will not provide a faithful de- 
scription of the monojet bounds if the limit on A is much below the pT cut. In the presence 
of a light mediator rj the relic density depends upon the both the mass tbdm of the DM 
state X and the mass of the mediator tti^. The annihilation of the symmetric component is 
enhanced in the mass range moM ^ ?7i^, as the annihilation can proceed by the t-channel 
processes XX — )• rjr]. Moreover, when tti^ ~ 2mDM the annihilation via the s-channel process 
XX — 7- ?7 — 7- SM is resonantly enhanced. Consequently, the case of GeV scale mediators must 
be carefully analysed and it is not sufficient to simply adjust the monojet limits as is often 
the approach taken in the previous literature. 

In the remainder of this section we use the package micrOMEGAs [61] to study two well 
motivated examples and demonstrate the full effect of light mediators on relic density cal- 
culations and experimental bounds. To relate the features of the micrOMEGAs results to the 
underlying physics, we study in the Appendix, a semi-analytic approach, following [84]. 
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Figure 4. Constraints on a light scalar mediator coupling fermion DM to the visible sector via mixing 
with a SM-like Higgs with mass m/i = 125 GeV. The black curve shows the minimum DM-mediator 
coupling Ax required to efhciently annihilate the symmetric component of the DM. The red curve 
indicates the current combined spin independent direct detection bounds. Constraints from monojet 
searches are negligible. The allowed parameter space is indicated by the shaded region. For mediators 
heavier than a few GeV only the resonance region survives the tension between direct detection liinits 
and the requirement for efficient annihilation of the symmetric component. 



4.1 The light scalar-Higgs portal 

To illustrate this general formalism we shall consider a particular example of a scalar state 
which mediates interactions between fermion DM il) and SM quarks via mixing with the SM 
(or a SM-like) Higgs. In minimal models of DM it is expected that the mediator connecting 
the hidden and visible sectors should be a singlet under the SM gauge groups. In this case, 
mixing with the Higgs after electroweak symmetry breaking provides the dominant manner 
by which the mediator state couples to the SM quarks. After electroweak symmetry breaking 
the scalar mediator and the Higgs will mix giving physical mass eigenstates [j] and K) with 
interactions of r] with the SM quarks induced by this mixing [85, 86], which we parameterise 



£ D Xxrj^'fp + '^{OX'yq)m'l, 



(4.1) 



where yq are the standard model Yukawa couplings, 9 



mu 



is the mixing parameter and 



the product of these may be dressed by an additional coupling A'. We use micrOMEGAs to 
calculate the required \x necessary to obtain the asymmetry dominated relic density, as a 
function of DM mass w-dm- See the Appendix for a semi-analytic treatment of ADM relic 
density in the presence of s-channel resonances il^ijj — )■ r; — t- SM and t-channel threshold 
effects ifjil) — >• rjrj. The results are shown in Fig. 4 for a scalar mediator with with parameters^ 
{rn,^, A'} = {10 GeV, 1}, {50 GeV, 0.1} and are accompanied by current direct detection 
constraints. For mediators heavier than a few GeV only the resonance region survives the 
tension between theoretical and experimental requirements. Moreover, because of the strength 



®LEP searches [87] for mixed Singlet-Higgs states constrain \'6 to be small and influence are choice of 
values for A', however varying this parameter will only lead to rescaling of the results. 
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Figure 5. Constraints on a pseudoscalar mediator coupling fermion DM to the visible sector. The 
solid curve shows the required hidden sector coupling in order to efficiently annihilate the symmetric 
component. The dashed curve visible near 5 GeV (the resonance region) indicates the LHC monojet 
limits, as can be seen these do not present significant constraints on the model, this is in contrast to 
the corresponding contact operator O^ studied in Section 3. Moreover, as the scattering cross-section 
is suppressed by g^ there are essentially no limits from direct detection. The allowed parameter space 
is indicated by the shaded region. 



of the spin-independent direct searches even lighter mediators are constrained to the point 
where only a small part of the parameter space survives. For scalar mediators with mass < few 



GeV, for twdm 



m 



r], 



DM annihilation into rj pairs (which later decay to SM states) is allowed 



by direct detection constraints and presents a possible mechanism. However, mediators with 
"^r; ^ 10 GeV or DM with moM ^ 5 GeV will lead to additional limits, most prominently 
searches for quarkonium decays to invisible states [74]. If the experimental hints of a SM-like 
Higgs boson at approximately 125 GeV are confirmed, then this portal interaction will be 
constrained by subsequent precision measurements of Higgs decay rates [88-90] . 

4.2 The light pseudoscalar portal 

It was seen in Section 3 that heavy mediators are greatly disfavoured in natural implemen- 
tations of ADM. Pseudoscalars are an ideal candidate for light mediator states and present 
a particularly interesting case as the constraints from direct detection (suppressed by g^) are 
negligible. We parametrise the portal interaction connecting fermion DM ip with SM quarks 
via a pseudoscalar mediator as follows 



£ D i\xaip-f^ip + ^ i{0X'yg)aq^^fq, 



(4.2) 



where, as previously, yq are the standard model Yukawa couplings and 9 is the mixing pa- 
rameter. Such a scenario could occur via mixing with a CP-odd state of an extended Higgs 
sector (e.g. the A° of the MSSM). 

In Fig. 5 we display (black solid curve) the minimum Xx required for efficient annihilation 
of the symmetric component, as a function of itidm, for a portal interaction due to a 10 GeV 
pseudoscalar mediator a for 9\' = 0.1, 0.01 as the black solid curve. By fixing the coupling 
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of the mediator to SM quarks, the size of the DM couphng Xx to the mediator required 
for efficient annihilation of the symmetric component is simply a function of the DM mass 
and the mediator mass. This allows us to present the effects of varying the parameters in a 
clear manner, which can be readily compared with the previous example of the light scalar- 
Higgs portal. The leading constraints on the light pseudoscalar portal come from the LHC 
monojet searches. We use CalcHEP to determine the LHC monojet bounds following the same 
procedure as in Section 3 and these limits are presented as the red dashed curve in Fig. 5. 

We note that both the monojet limits and the DM relic density are very sensitive to the 
DM mass, with both curves featuring resonances and mass threshold effects. In contrast to 
the effective operator analysis, where ADM was excluded up to 100 GeV by collider searches, 
the monojets limits for light mediators nia ^100 GeV are weakened to the extent that for 
natural couplings they present no appreciable constraint. Indeed the monojet limits are only 
visible in the plots of Fig. 5 in the resonance region nia ~ 2mDM- 

For iTia 7^ 2?7i£)M the DM annihilation is via a virtual mediator and this leads to a large 
suppression in the cross-section, thus for natural values of the parameter the mediator and 
DM masses must typically be arranged such that the annihilation of the symmetric component 
proceeds via resonant s-channel processes or t-channel pair production of the mediator states. 
In particular, if the DM is heavier than the mediator then there will generally be unconstrained 
parameter space in which efficient annihilation of the symmetric component to the visible 
sector can be achieved. Notably, light pseudoscalars states are very well motivated, as they 
typically occur as pseudo-Nambu-Goldstone Bosons. 

5 Concluding remarks 

If one begins from the well motivated assumptions that ADM with mass moM ~ '^p anni- 
hilates dominantly to SM quarks, and the flavour structure of such DM-quark couplings is 
minimal, then we have shown (Figures 1-3) that the current experimental limits allow one to 
make decisive statements regarding the nature of the hidden sector. In Fig. 1 it can be seen 
that contact operators with spin-independent non-suppressed direct detection cross-sections 
are excluded for A <1 TeV with universal couplings and up to 200 GeV in the case of ruq 
dependent couplings. The one exception is the operator \(t)^4>qq for m-DM ;$ 1 GeV. We 
expect that this small window will be explored by ATLAS and CMS in their 2012 run, where 
> 10fb~ will be collected by each experiment. Contact operators with spin-dependent non- 
suppressed direct detection cross-sections are displayed in Fig. 2 and it is seen that these are 
excluded by monojet searches up to A ~ 300 GeV. Finally, the suppressed direct detection 
cross-sections displayed in Fig. 3 are excluded for A <100 GeV with the exception of the 
pseudoscalar operator '^ip'f^ipq'j^q which is ruled out in the range A <10 GeV. 

Collating these results we conclude that most contact operators are very strongly con- 
strained and, moreover, that all of the operators are essentially excluded in the interesting, 
motivated ADM mass region 1 GeV < "idm ^ 10 GeV. However, the effective operator 
description is no longer valid for mediators lighter than <100 GeV, at which point resonance 
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and threshold effects become important in determining the exclusion limits. As can be seen 
in Figures 4 and 5, if the mediator states are light then successful models of ADM can be 
constructed. Hence, because of the strong constraints on contact operators, it is required 
that the mediator connecting the visible and hidden sectors be relatively light in models of 
ADM conforming with our natural assumptions in order to efficiently annihilate the symmet- 
ric component whilst remaining consistent with search constraints. Moreover, it can be seen 
in Figures 4 and 5 that for mediators with mass ~ 10 — 100 GeV there remain strong limits 
in many cases and to satisfy these constraints it is often required that the DM and mediator 
masses are arranged such that the DM annihilation is resonantly enhanced. 

Some caveats to the conclusion presented above are in order. Most prominently, the 
DM may couple dominantly to leptons. Although such a coupling would go against our 
expectations from models of flavour and from grand unified theories, if the DM interacts 
more strongly with leptons compared to quarks, this would greatly relax the monojet bounds 
and the constraints from DM direct detection [91]. The dominant bounds would now come 
from monophoton experiments and DM experiments which do not veto on electron recoil, such 
as DAMA/LIBRA [92]. Alternatively one may attribute a non- minimal flavour structure to 
the interaction with quarks, for instance, the DM could couple in an isospin violating fashion 
[93]. There exist studies of the collider limits in both the leptophilic [45] and isospin violating 
[47] scenarios with reference to conventional models of DM. Whilst possible in the context of 
ADM models these variant theories require further model building and a dedicated study. 

The introduction of new light states is perhaps the most interesting and well motivated 
manner to avoid the strong constraints presented here. If the DM belongs to an extended 
hidden sector, containing lighter states into which the symmetric component can annihilate, 
then the limits from both direct search constraints and cosmological requirements are greatly 
ameliorated. Indeed, entire hidden sectors are well motivated in many models of physics 
beyond the standard model. In particular, it has been argued that the topological complexity 
of generic string compactifications leads to the expectation of many hidden sectors sequestered 
from the SM states [94-97] . Studies of hidden sector cosmology and collider phenomenology 
exist in the literature [98-101] and this is also related to the hidden valley scenario [102]. 

The findings of this paper show that experimental searches strongly motivate the existence 
of additional hidden light degrees of freedom in models of ADM, either to mediate interactions 
between hidden and visible sector or into which the DM can annihilate. If the mediator state 
is made sufficiently light then portal interactions can provide mechanisms for ADM that both 
provide efficient annihilation of the symmetric component and satisfy limits from monojet and 
direct detection searches. However, depending on the mass and nature of the mediator there 
are further model dependent constraints on the couplings of these states, such as invisible 
quarkonium decay, beam dump experiments, invisible Z decay width, electroweak precision 
measurements, muon g — 2, atomic parity violation, and also astrophysical processes. The 
cosmology, phenomenology and model building opportunities presented by extended hidden 
sectors of ADM is the focus of our forthcoming companion paper [55]. 
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A Semi-analytic analysis of light mediators 

We shall make some remarks regarding the source of the features exhibited in the relic density 
calculations with light mediators using a semi-analytic analysis. The introduction of light 
mediators leads to resonance and mass threshold effects and the dominant process which 
determines the relic density becomes strongly dependent upon the DM mass. Note also that 
many models, in particular supersymmetric models, may lead to co-annihilation effects which 
must be accounted for but this is beyond the scope of this short appendix. In correctly 
specifying the annihilation integral J for light mediators, we shall follow the work of Griest 
and Seckel [84]. 

With a light mediator i] the standard expansion of {av) used for contact operators, eqn. 
(2.4), is not valid over the whole mass range of the DM, X. There are a number of notable 
ranges of the DM mass in which the annihilation of the symmetric component is enhanced. 
In the DM mass range moM ^ "^^q the annihilation can proceed by the t-channel process 
XX — )• r/77. The second region of interest is when rriri ~ 2m,DM) in which case annihilation via 
s-channel XX — )• 77 — )• gg is resonantly enhanced. 

To model the effects of resonant production we factor out the pole factor P{v^) from the 
cross-section 



P(v^) 



2\ ^ / n \ 2 



(2mDM/m^) \ ,f^ 



1-T 



(A.l) 



Whilst more sophisticated approximations are discussed in [84], this will suffice for our pur- 
poses. Following this we can write the thermally averaged cross-section in the mass region 

"^DM ^ m-r] as follows 

i^^)muM<m., ^ [as + ^) X P(0), (A.2) 

where a^ and bg are the coefficients from the velocity expansion of the cross-section correspond- 
ing to s-channel XX — )• r/ — )• qq. It is clear from inspection of eqn (A.l) that the cross-section 
will be peaked at rriDM ~ m^/2. In which case the annihilation integral, appearing in eqn. 
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(2.5) may be expressed as follows 

'as + f^]x P(0) 



'^rnuM<mr, ~ ^ ■ (^-3) 

In the mass region rrir] > m-DM the states rj can be pair produced and present the dominant 
annihilation channel. For DM states with mass ttidm ^ "ir; the final state velocity vj becomes 
important and can not be approximated as unity in the velocity expansion of the cross-section 

av ^ {at + btv'^)vf. (A. 4) 

The quantity Vf can be expressed as follows 



We may make a velocity expansion of Vf 




ml ( mlv'^ \ 

-^^\^ + ^r^ 2Y • (A-6) 

Note that for m^ ^ "idm the vj ~ \ and we recover the standard velocity expansion for 
the cross-section. Furthermore, observe that the velocity expansion for Vf breaks down for 
'TT'DM ~ t^'q- The annihilation integral for ttidm > "ir; and away from the threshold region 
"T-DM ~ T^rii may be approximated thus 

7- ^ ~ V ^^ / /A 7^ 

where at and ht are the coefficients from the velocity expansion of the cross-section corre- 
sponding to t-channel XX — )• 7]%]. The factors la^b account for the final state velocities of the 
produced states and are given by: 

4=^r<^d., and 4=?£lr<5^d.. (A,8) 

'^i Jxp ^ '^t Jxp OX 

The explicit form of the factors Ia,b given in eqn. (A. 8) depend on Vf. Away from the threshold 
region the approximation in eqn. (A. 6) is valid, however near the threshold one must use the 
exact form for Vf. Moreover, since the interacting states have non-zero velocities, processes 
which would be forbidden at zero velocity may be accessible and lead to the dominant decay 
channels. We construct the total annihilation integral as a piecewise function: 



J 



'Jrnj^M<m^ "^DM < 0.97 X m^ 

-^Threshold 0.97 X nir, < moM < xl.OSm^ , (A. 9) 
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where Jrhreshoid is identical in form to JmuM>m but the la^b factors are evaluated using the 
exact form for Vf. The factors /„,& ^ire given explicitly in [84] for all cases. 

To calculate the relic density with light mediators we must use the corrected annihilation 
integral in eqn. (2.1): 

^DM = Z f ;— r, -fn 



1 - exp [-riDuJ^] ' °^ exp [rj^uJuj] - 1 ' 

Then, as previously, we require that the asymmetric component alone accounts for the relic 
density by the following requirement on the symmetric component: 



J_ 17dm fe^ /GeV\ 



>Sym < TT^X 



At this stage, for a given model, one can numerically evaluate the analytic forms above. 
This semi-analytic computation gives comparable results to those obtained from micrOMEGAs. 
Consequently, the features appearing in the light mediator analysis of Section 4 can be well 
understood by comparison to the semi-analytic approach. 
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